Ultracold dipolar gases - a challenge for experiments and theory 
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We present a review of recent results concerning the physics of ultracold trapped dipolar gases. In 
particular, we discuss the Bose-Einstein condensation for dipolar Bose gases and the BCS transition 
for dipolar Fermi gases. In both cases we stress the dominant role of the trap geometry in determining 
the properties of the system. We present also results concerning bosonic dipolar gases in optical 
lattices and the possibility of obtaining variety of different quantum phases in such case. Finally, 
we analyze various possible routes towards achieving ultracold dipolar gases. 



I. INTRODUCTION 

In the last twenty years at least three Nobel Prizes were 
awarded for studies of the phenomena of condensation, 
superfluidity and superconductivity 0, whereas several 
others were also related to these subjects [0. The three 
Nobel Prizes in atomic, molecular and optical (AMO) 
physics within the same period mark the path from trap- 
ping and cooling of ions, through laser cooling of atoms, 
towards Bose-Einstein condensation (BEC) of trapped 
atomic gases . The observation of BEC || has bursted 
a new interdisciplinary area of modern AMO and con- 
densed matter physics: the physics of ultracold weakly 
interacting trapped quantum gases . So far, most of the 
experimental results in this area can be very accurately 
modeled by mean- field methods and its extensions, based 
on the Gross-Pitaevskii (GP) equation and Bogoliubov- 
de Gennes equations for bosonic gases ||, and on the 
BCS theory for fermionic ones ||. 

This new research area presents several novel aspects, 
in particular the finite character and inhomogeneity of 
the considered systems, and perhaps more important, the 
unprecedented possibility of control and manipulation of 
the system properties, which allows to study situations 
that were not yet encountered in condensed matter or low 
temperature physics. The number of parameters which 
can be successfully controlled is large, e.g. temperature, 
number of atoms, trap potentials, interatomic interac- 
tions, etc. In addition, the internal level structure of the 
atoms can also be employed, e.g. to manipulate BECs 
using Raman-Bragg techniques M, or to analyze multi- 
component condensates H 



*This paper is based on the lecture given by M. Lewenstein 
at the Nobel Symposium "Coherence and Condensation in 
Quantum Systems", Gothesburg, 4-7.12.2001. 



It is also worth mentioning, that the experimental tech- 
niques have progressed recently to a stage at which mean 
field methods cease to provide a proper physical picture. 
In this sense, recent experiments at JILA S in which 
the scattering length can be modified at will by using 
Feshbach resonances, allow to realize systems with very 
large scattering length, in which the mean-field picture is 
no more applicable. Similarly, the realization of a Bosc- 
condensed metastable Helium gas |h|, with a potential 
to study higher order correlation functions of the sys- 
tem, demands an analysis beyond the mean field theory. 
The recent demonstration of the Mott insulator to su- 
pcrfluid phase transition with atoms in an optical lattice 
]TT| , predicted in Ref. , belongs to the same category, 
but at the same time opens a new research area of AMO 
physics: the physics of strongly correlated quantum gases. 
The system in question allows for an easy and accurate 
control and manipulation, and thus provides a new and 
particularly promising test ground for theories of quan- 
tum phase transitions Jl3| , which have traditionally dealt 
with condensed-matter rather than atomic systems. 

In recent years considerable interest has been devoted 
to another aspect of the internal structure of the parti- 
cles forming an ultracold gas, namely their dipole mo- 
ment. If such dipole moment is sufficiently large, the re- 
sulting dipole-dipole forces may influence, or even com- 
pletely change the properties of BEC in bosonic gases, 
the conditions for BCS transition in fermionic gases, or 
the phase diagram for quantum phase transitions for ul- 
tracold dipolar gases in optical lattices. In this paper we 
present a review of our recent results on ultracold dipolar 
gases. We consider here only the case of polarized dipolar 
gases when all dipoles are oriented in the same direction, 
and discuss the effects of the dipole-dipole interactions 
and their interplay with short range (Van der Waals) in- 
teractions. We do not consider several other important 
aspects of dipole-dipole interactions, e.g. the role of the 
dipole-dipole interactions in the spontaneous polarization 
of spinor condensates in optical lattices [Q , or the self- 
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bound structures in the field of a traveling wave . As 
dipolar interactions can be quite strong (relative to the 
typical low-energy collisions characterized by the s-wave 
scattering length), the dipolar particles are considered as 
promising candidates for the implementation of fast and 
robust quantum-computing schemes Jl(],[l7j ■ 

Our review is organized as follows. Section || is de- 
voted to bosonic dipolar gases, the prospects for their 
condensation, the ground-state properties of the BEC, 
and its elementary excitations in a trap. In sections 
III, an analogous analysis is applied to ultracold dipo- 
lar Fermi gases. In particular, the possibility of the BCS 
transition is investigated. Section IV briefly considers 
the issue of quantum phase transitions in a dipolar Bose 
gas placed in an optical lattice. In section [v| we discuss 
various possible physical systems which could lead to an 
ultracold dipolar gas. Finally, we conclude in section jv|. 



II. ULTRACOLD DIPOLAR BOSE GASES 

In this section, we consider a system of N bosonic par- 
ticles possessing a dipole moment, and confined in a har- 
monic trap of cylindrical symmetry. We constraint our 
analysis to the case in which all dipoles are assumed to 
be oriented along the symmetry axis of the trap. Addi- 
tionally, we assume that the particles interact via dipole- 
dipole forces, and that these forces either play a dominant 
role, or at least compete with the short-range forces. 

Why are dipolar gases interesting? The answer is sim- 
ple: because they are dipolar. The dipole-dipolc inter- 
action potential between two dipolar particles is given 
by 

Vrf(R) = (d 2 /R 3 )(l - 3cos 2 6>), 

where d characterizes the dipole moment, R is the vector 
between the dipoles (R = |R| being its length), and 9 
the angle between R and the dipole orientation (Fig. 1). 
The potential Vd(R) has two important properties: it is 
anisotropic, and is of long-range character. As we discuss 
below, these properties have important consequences. 

At low temperatures one expects the dipolar Bose 
gases to condense. One expects also that the condensate 
properties will dramatically depend on the geometry of 
the trap. In cigar-shape traps along the dipole direc- 
tion the interactions will be mainly attractive, and the 
condensate will be unstable, similar to the case of a gas 
with attractive short-range interactions (negative s-wave 
scattering length) Q . Conversely, in pancake traps the 
interactions will be mainly repulsive, and the gas might 
become stable. Therefore, the dipolar gases offer the un- 
precedented possibility of modifying the atom-atom in- 
teractions by turning an "easy knob" , namely the trap 
geometry, which is relatively easy to control and modify 
exp er iment ally. 
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FIG. 1. Anisotropy of dipole-dipole interactions. 

Several groups have recently studied the physics of ul- 
tracold dipolar bosonic gases. L. You and his collabora- 
tors |h] |2l| have considered the question of the validity 
of the mean-field approach in this case. This is a non 
trivial question, since contrary to the case of short-range 
potentials, where the s-wave scattering always becomes 
dominant at very low temperatures, the scattering ampli- 
tude for dipole-dipole interactions has non vanishing con- 
tributions of all partial waves at low energies. Moreover, 
the interplay play between the dipole interactions and 
the short range forces might lead to shape resonances in 
the atom-atom scattering. Close to such resonances the 
effective scattering length diverges, and the mean field 
approach cannot be used. Away from shape resonances, 
however, the Gross-Pitaevskii (CP) equation, including 
the non-local dipole-dipole interactions, provides a good 
description of the condensate. Both L. You and his co- 
workers p9| , Coral et ol. and Mackie et al. [^3| have 
investigated the properties of the ground state of these 
systems, and the interplay between dipolar forces and 
the short-range interactions. We have concentrated in 
our studies [Q on the case of dominant dipole-dipole in- 
teractions. This case is discussed in the next subsection. 



A. Ground state properties of dipolar Bose gases 

Similarly to |l9|,|2| , we describe the dynamics of the 
condensate wave function tp(r,t) by using the time- 
dependent GPE: 
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Here ?/>(r, t) is normalized to the total number of conden- 
sate particles N . The third term in the rhs is the mean- 
field corresponding to the short-range forces, whereas the 
last term is the mean field of the dipole-dipole interac- 
tion. In the following, we omit the term g\t/}(r, t)\ 2 t(j(r, t), 
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assuming that the interparticle interaction is dominated 
by the dipole-dipole forces (d 2 3> \g\ = 4:irh 2 \a\/m, where 
a is the s-wave scattering length, and m is the parti- 
cle mass), and that the system is away from shape res- 
onances of Vd(R). The ground-state properties are gov- 
erned by the stationary GPE, in which the lhs of Eq. 
(|l|) is replaced by fiip(r), where the chemical potential /j, 
corresponds to the minimal energy solution. 

There are two important parameters that can be easily 
controlled in experiments: the collective dipole strength 
Nd 2 , and the trap aspect ratio I — {up/ojz) 1 / 2 — a z /a p , 
where the characteristic harmonic oscillator length ai = 
yjfi/muii. The first of these parameters can be rewrit- 
ten into a dimensionless form a — (Afd 2 /a^ ax )/?iw max , 
where a max = (h/2Muj rn [ ri ) 1 ^ 2 , with uj m i n the minimum 
trap frequency. Therefore, a represents the ratio of the 
dipole-dipole interaction energy at the characteristic har- 
monic oscillator distance to the characteristic harmonic 
oscillator energy, to rn i n . Another quantity, which criti- 
cally characterized the state of the system is the mean 
dipole-dipole interaction energy per particle given by the 
expression V = (1/jV) J_V d (r - r')V'o( r )'0o( r ') c!rrfr '- 

The results of Ref. J24J] can be summarized as follows. 
For cigar shaped traps with I > 1 the mean-field dipole- 
dipole interaction is always attractive, and the gas be- 
comes always unstable if the number of particles N ex- 
ceeds a critical value N c , which depends only on the trap 
aspect ratio I. The quantity \V\ increases with N and 
the shape of the cloud changes. In spherical traps the 
cloud becomes more elongated in the axial direction and 
near N = N c the shape of the cloud is close to Gaussian 
with the cloud aspect ratio L = L z /L p ~ 2.1. In cigar- 
shaped traps (I 3> 1) especially interesting is the regime 
where %uj z <C \V\ <C hu p . In this case the radial shape 
of the cloud remains the same Gaussian as in a non- 
interacting gas, but the axial behavior of the condensate 
will be governed by the dipole-dipole interaction which 
acquires a quasi 1-dimensional (ID) character. Thus, one 
has a (quasi) ID gas with attractive interparticle interac- 
tions, i.e. a stable (bright) soliton-likc condensate, where 
attractive forces are compensated by the kinetic energy 
p5j. With increasing N, L z decreases. Near N = N c , 
where |V| is close to Tuj p , the axial shape of the cloud also 
becomes Gaussian and the aspect ratio takes the value 
L « 3.0. 

The situation is quite different for pancake traps (I < 
1); in particular there exists a critical aspect ratio Z* ~ 
0.41, which splits the pancake traps into two different 
categories: soft pancake traps (Z* < I < 1) and hard 
pancake traps (I < £»). For soft pancake traps the dipole- 
dipole interaction energy is positive for a small number of 
particles and increases with N. The quantity V reaches 
a maximum, and a further increase in N reduces V and 
makes the cloud less pancake. For a critical number of 
particles N — N c the BEC becomes unstable. 

We have found generally that the dipolar condensate 



is unstable and collapses when N > N c for V < with 
|V| > huj p . However, for hard pancake traps with I < I* 
the condensate is stable at any N, because V always re- 
mains positive. For small ./V the shape of the cloud is 
Gaussian in all directions. With increasing N, the quan- 
tity V increases and the cloud first becomes Thomas- 
Fermi in the radial direction and then, for a very large 
N, also axially. The ratio of the axial to radial size of the 
cloud, L — L z / L p , continuously decreases with increas- 
ing number of particles and reaches a limiting value at 
N — > oo. In this respect, for a very large N we have a 
pancake Thomas-Fermi condensate. 

It is worth stressing that many of the above results, 
calculated from a direct numerical simulation of the time- 
independent GPE, were also analytically obtained with 
the help of a variational method already used in the con- 
text of short-range interacting BEC (c.f. p|Q). 
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FIG. 2. Graphical representation of oscillations modes of 
the condensate. 



B. Excitations in dipolar Bose gases 

In this section we analyze the elementary excitations 
in a trapped dipolar gas. One possible method of analy- 
sis of these excitations could be to solve the correspond- 
ing Bogoliubov-de Gennes equations pj, which become 
non-local due to the presence of dipole-dipole interac- 
tions . Another approach, suggested in the context of 
short-range interacting condensates by Ruprecht et at 
p8| , is to use the time-dependent GPE, and numerically 
study the spectrum of the small perturbations around 
the ground-state solution. For the lowest lying excitation 
modes, such those schematically sketched in Fig. H, it is 
also possible to obtain analytic results using the time de- 
pendent variational method of Ref. [p9| , which provides 
very accurate results when a is not to close to its crit- 
ical value cr c , which is defined as the dipole strength at 
which the BEC becomes unstable. L. You and his col- 
laborators [^0],^l| have recently employed this method to 
analyze various properties of the excitations in dipolar 
Bose gases. We have recently used the numerical 
method of Ref. |2j| and the variational approach of Ref. 
p9[ , in order to answer two fundamental questions con- 
cerning elementary excitations, namely how is the qual- 
itative nature of the instability, and how the effects of 
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the dipole-dipole interactions can be observable from the 
excitation spectrum. 

The question of the nature of instability is twofold. 
First, the determination of which of the modes becomes 
unstable when the number of atoms reaches the critical 
value. Second, the analysis of the behavior of the cor- 
responding mode frequency close to the criticality. This 
analysis has been first performed for the case of a Bose 
gas with attractive short-range interactions. Bergeman 
pl[ observed from his numerical results that, as the ra- 
tio of the nonlinear interaction energy to the trap fre- 
quency, 7, approaches the critical value, the frequency of 
the "breathing" mode 2 tends to zero and merges with 
the frequency of the Goldstone mode corresponding to 
the overall phase of the condensate. The parameter 7 is 
defined as 
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where cto is the geometric mean width of the ground state 
in an ellipsoidal harmonic trap. In the interesting region 
of a < 0, the critical value of 7 is j c = —0.54. Above 
criticality, the "breathing" becomes unstable and attains 
complex frequency. Singh and Rokhsar analyzed this in- 
stability using self-similar solutions to describe the 
modes, or equivalently the variational approach of the 
previous subsection. They have shown that close to crit- 
icality the frequency of the mode 2 vanishes as I7 — 7 C | 1 ^ 4 . 

In the case of dipolar gases with dominant dipole in- 
teractions the situation is completely different and much 
more complex. Only for aspect ratios far above the crit- 
icality, I 3> I* (I > 1.29) the situation reminds that of a 
gas with attractive short-range interactions. The mode 
corresponding to the lowest frequency is the "breathing" 
mode 2. This mode becomes unstable when the param- 
eter (j — ► cr c . The scaling behavior of its frequency u>2 
can be analyzed employing the variational approach of 
the previous section, and the analytic (although approx- 
imate) expressions for U2 from Ref . [2(J . We find that ui2 
goes to zero as (<r c — a) 13 , with (3 ~ 1/4. 

For intermediate values of I > I* (0.75 < I < 1.29), 
the exponent j3 is still close to 1/4, but the "breathing" 
and quadrupole modes mix. For a far below a c the mode 
corresponding to the lowest frequency is the "breathing" 
mode 2. As we approach the critical value of a the char- 
acter of the lowest frequency mode changes and becomes 
quadrupole- like. 

For I close to I* (I < 0.75) the situation changes and 
the mode corresponding to the lowest frequency is the 
quadrupole mode 3. Now, it is its frequency u>3 which 
tends to zero as the parameter a approaches the critical 
value. For / not too close to I* the exponent (3 is still 
close to 1/4. Completely new effects arise due to the 
existence of the previously discussed critical aspect ratio 
I* ~ 0.41. This is illustrated in Fig. |^. As one approaches 
/*, the exponent /3 departs from 1/4 towards a greater 



value. This crossover is explained in Fig. [| For I slightly 
below I*, the frequency of the quadrupole mode W3 has a 
quadratic minimum close a = a c (l*). Exactly at I = I*, 
0)3 goes thus to zero as (a c — a) 2 , i.e. (3 = 2. 




FIG. 3. Frequency 013 of the lowest mode 3 as a function of 
a for I just below (solid line), and just above (dashed line) /*. 
For I < I* , the frequency UI3 has a quadratic minimum close 
to a = a c (l*). 




trap aspect ratio co_/(0 

FIG. 4. Difference of excitation frequency for modes 1 
(dashed line), 2 (dotted line) and 3 (solid line) between the 
cases of purely contact and mixed (contact and dipole-dipole) 
interactions as a function of the trap aspect ratio for 10000 
atoms and a — as a (data from variational analysis). 

Apart from the question of the nature of instability, 
in Ref. |5(| we have considered also the situation in 
which the particles interact via both contact and dipole- 
dipole forces. In this case the strength for both types 
of forces is comparable and neither of them can be ne- 
glected. As a specific example we have considered the 
case of 52 Cr, which has drawn some experimental in- 
terest |33pB6[. Chromium has a large magnetic dipole 
moment of 6/is (Bohr magnetons), but its s -wave scat- 
tering length is still unknown. We have considered the 
trap frequencies of an experiment underway at the Uni- 
versity of Stuttgart |Q: uj z = 2n40 Hz and uj p = 27r485 
Hz. Our results suggest clearly that the effects of dipole- 
dipole interaction would be accessible to experimental de- 
tection if condensation of chromium or europium p8j is 
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achieved, and if their scattering lengths are not anoma- 
lously large. Dipole-dipole effects can be enhanced by 
reducing the s-wave scattering by employing Feshbach 
resonances [|9||4(J. It is thus very important to analyze 
possible parameters (scattering length, number of atoms, 
trap aspect ratio) that would maximize the predicted fre- 
quency shifts. A typical results of our analysis is shown 
in Fig. |. 



III. DIPOLAR FERMI GASES 

The recent success in observing the quantum degener- 
acy in ultra-cold atomic Fermi gases J|l]-Q stimulates 
a search for gaseous Fermi systems where combined ef- 
fects of Fermi statistics and interparticle interactions re- 
sult in a non-trivial physical behavior. Due to the Pauli 
principle, to observe these effects in the case of a short- 
range Van der Waals interparticle interaction, a simulta- 
neous trapping of at least two different fermionic species 
is needed, with a rather severe constraint on their relative 
concentrations. The situation is different for Fermi gases 
of dipolar particles. In the ultra-cold limit, the dipole- 
dipole scattering amplitude is energy independent for any 
orbital angular momenta L > 0. This observation follows 
from the results of p5| , j46| . It has also been recently dis- 
cussed with important physical insight and consequences 
by You and coworkers [}47 48 1. This opens prospects to 



observe the effects of interparticle interaction in a single- 
component Fermi gas, where only scattering with odd 
orbital momenta (negligible in the case of Van der Waals 
interactions) is present. These prospects are especially 
interesting as in single-component fermionic gases the 
Pauli exclusion principle provides a strong suppression 
of inelastic collisional rates (see 0]). Hence, one can 
think of achieving significantly higher densities than in 
Bose gases. 



A. Physical system 

We consider a single component gas of fermionic dipo- 
lar particles in a harmonic trap, with dipole moments 
oriented in the same direction (z -axis). Then, due to 
the Pauli principle, the contribution of the short-range 
part of the interparticle interaction can be neglected, and, 
therefore, the dominant interaction between particles is 
the dipole-dipole one. As a result, the Hamiltonian of 
the system reads 

H = J drft(r) |-|^A + V trap (r) - M | ^(r) (3) 
+ i f dvdr'^(r)^(v)V d (r-r')^(r')^(r'), (4) 



where Vd(r) = (d 2 /r 3 )(l — 3 cos 2 (6' r )) is the dipole-dipole 
interaction with 9 r being the angle between the interpar- 
ticle distance r and the z-axis, n is the chemical potential, 
and Vt rap (r) = m[oj 2 z 2 +uj 2 (x 2 +y 2 )}/2 is the trapping po- 
tential; -0(r) and ip^{r) denote annihilation and creation 
operators for fermions with canonical anticommutation 
relations. 

We assume that the gas is in the regime of quan- 
tum degeneracy, i.e. the temperature T is much smaller 
than the chemical potential /i (or the Fermi energy Ef), 
T <C [i = £f = (fi. 2 /2m)(67r 2 no) 2 / 3 , where no is the max- 
imal gas density (in the center of the trap) , and the inter- 
particle interaction is weak. The latter means that the 
mean dipole-dipole interaction energy per particle nd 2 is 
much less than the Fermi energy, nd 2 <C Ep- 



B. Why are dipolar Fermi gases interesting? 

Which features make this system attractive and chal- 
lenging for both theorists and experimentalists? The sys- 
tem under consideration is unique in that the physics 
is dominated by the long-range dipole-dipole forces (the 
short-range part of the interparticle interaction, as al- 
ready mentioned, can be neglected due to the fermionic 
nature of the particles). The dipole-dipole interaction is 
anisotropic and, as a result, it mixes scattering channels 
with different angular momenta I. Additionally, the in- 
teraction is partially repulsive (when two particles are 
side by side to each other) and partially attractive (when 
they are on top of each other). This means that, similarly 
as for dipolar Bose gases, the properties of the system de- 
pend on the trap geometry (see Ref. |4|]). Indeed, in a 
pancake trap, most of the time the particles are side by 
side to each other, and the average interparticle interac- 
tion is repulsive. In the opposite case of a cigar shaped 
trap, the particles are mostly on top of each other and, as 
a result, the average interaction is attractive. Therefore, 
one can expect a pronounced dependence of the system 
properties on the aspect ratio I = ^/lo p /uj z of the trap. 

The fact that the dipole-dipole interparticle interaction 
is partially attractive opens the possibility for the BCS 
pairing at sufficiently low temperatures. For example, 
the matrix element (L = 1, M = 0\V d \ L = 1, M = 0) of 
the dipole-dipole interaction between states with angular 
momentum L = 1 and its projection to the z-axis M = 0, 
is negative: 

47T 

(L = 1, M = |Vd| L = l, M = 0) = d 2 < 0. 
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This signals about the possible BCS transition 50 52 . 
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C. The BCS pairing 

The BCS pairing transition corresponds to the for- 
mation of Cooper pairs, the correlated states of two 
fermions. The corresponding order parameter A that 
appears below the transition temperature T c , describes 
the coherent motion of the Cooper pairs: 

A(r 1; r 2 ) oc (^(n)^(r 2 )) ^ 0. 

The order parameter vanishes for T — > T c and obeys 
the gap equation. For temperatures close to T c this equa- 
tion has the form (the so-called linearized gap equation) : 

A(ri, r 2 ) = -Vd(ri-r 2 ) / rfr 3 dr 4 ]i'(r ll r 2 ; r 3 , r 4 )A(r 3 , r 4 ), 

(5) 

where the kernel K is defined as 

tanh(^ 1 /2T) + tanh(£, 2 /2T) 



^(ri,r 2 ;r 3 ,r 4 ) = ^ 



<t> Vl (ri)^ 2 (r2)d ( r s)^ 2 (1-4) 

with £„ = e v — [i and 4> v (r) being the solutions of the 
single-particle Schrodinger equation 



- — A + F trap (r) \ 4> u (r) = e v (j) v {T). 



The critical temperature T c can be defined as the max- 
imum temperature for which the linearized gap equation 
(|J) has a nontrivial solution. It should be mentioned 
that this equation can only be used for finding the criti- 
cal temperature and the spatial dependence of the order 
parameter. To determine the absolute value of the order 
parameter, which is temperature dependent, one has to 
consider the terms nonlinear in A, omitted in Eq.(||). 



D. Spatially homogeneous dipolar Fermi gases 

In the spatially homogeneous case (Vt ra p(r) =0), the 
order parameter depends only on the relative coordinate 
r = ri — r 2 , A(ri,r 2 ) = A(ri— r 2 ), and the gap equation 
(||) with the leading non-linear term added, was discussed 
in Ref. f3^ ] , where it was found that the critical tempera- 
ture of the BCS transition in the spatially homogeneous 
Fermi gas of dipolar particles equals (including the so- 
called Gor'kov and Melik-Barkhudarov corrections) 

7T£p , . irh 

T c = l.Ue F exp(--— 2) = 1.44e F exp(-— , 

\2nd l 2\ad\PF 



(6) 



where 



is the gas density, ef 
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(h 2 /2m)(67r 2 n) 2 / 3 the Fermi energy, pp the Fermi mo- 
mentum, and ad = —2md 2 /7r 2 H 2 the effective scattering 



length. The latter has been introduced in order to make 
the exponent in the expression for T c , Eq.(||), look like the 
exponent in the expression for the critical temperature 
in the case of a two-component Fermi gas with the inter- 
component s-wave scattering length a<j. For the case of 
a single-component gas of fermionic ND3 molecules with 
the dipole moment d = 1.5D, the transition temperature 
T c is larger than 100 nK at densities n > 5 • 10 12 cm~ 3 . 

The properties of the superfluid dipolar Fermi gas are 
different from those of a two-component fermionic gas 
with s-wave pairing due to a short-range inter-component 
interaction. In the case of the s-wave paring, the order 
parameter is isotropic, whereas it is anisotropic in the 
superfluid dipolar Fermi gas. In the latter case, the value 
of the order parameter A at the Fermi surface p = pp , is 

A(p F , 9) = 2.5T cV /(T c - T)/T c ■ V2 sin(| cos(0)), 

where 9 is the angle between the momentum p and the 
z-axis. As a result, the gap in the spectrum of single- 
particle excitations, which appears below the transition 
temperature T c , is anisotropic. For example, excitations 
with momenta in the direction of the dipoles acquire the 
largest gap. In contrast to this, the eigenenergies of ex- 
citations with momenta perpendicular to the dipoles re- 
main unchanged. The properties of collective excitations 
are also expected to be dependent on the direction of 
their momenta. Therefore, the response of the dipolar 
superfluid Fermi gas to small external perturbations will 
have a pronounced anisotropic character. 

Another distinguishing feature of the superfluid dipo- 
lar Fermi gas is related to the temperature dependence of 
the specific heat. Well below the critical temperature the 
single-particle contribution to the specific heat is propor- 
tional to T 2 , rather than being exponentially small as in 
the case of the s-wave pairing. This follows from the fact 
that the energy e of single-particle excitations has a line 
of zeros on the Fermi surface: e(pf ) = for the angles 
at which A(pjr,n) = 0, i.e. for 9 = tt/2 and an arbi- 
trary azimuthal angle ip. As a consequence, the density 
of states in the vicinity of the Fermi energy is v{e) ~ e for 
e -C A . Therefore, at temperatures T -C A ~ T c , the 
temperature dependent part of the energy of the system 
is proportional to T 3 , and the specific heat is thus pro- 
portional to T 2 . This contribution is much larger than 
that of collective modes, which is oc T 3 and is dominant 
in the case of the s-wave pairing. 



E. BCS pairing in a harmonic trap 

We now discuss how the results of the previous Section 
change when the harmonic trapping potential is switched 
on. It turns out (see Ref. J53| for details), that the pres- 
ence of the trapping potential always results in the de- 
crease of the critical temperature of the superfluid tran- 
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sition. In the case where the trap frequencies u> z , oj p are 
much smaller than the critical temperature T c , Eq.(^|), 
of the transition in the spatially homogeneous gas with 
density n equal to the maximum gas density no in the 
trap, one obtains 



^trap rp 




where T c trap is the transition temperature in the trap, 
£(z) the Riemann zeta-function, and w = (w z uj 1 p ) v l' i . One 
can see from this expression that for sufficiently small 
ui/T c , the critical temperature in the trap T* rap is only 
slightly lower than T c . As follows from Eq.(Q), the opti- 
mal value of the trap aspect ratio I, that corresponds to 
maximal T c tra P at a given T c , is I ~ 1.38. 

Another interesting feature of the BCS transition in a 
trapped Fermi gas of dipolar particles is the existence of 
the critical value for lo z . This comes from the fact that 
the paired states in the trapped dipolar Fermi gas have 
different quantum number n z , and, therefore, their en- 
ergies are always different, at least by the amount oj z . 
When the difference becomes of the order of the order 
parameter A ~ T c which measures the strength of the 
paring correlations, the pairing is obviously impossible. 
As a result, the superfluid transition in the dipolar Fermi 
gas is possible only in traps with u z < uj zc , where the 
critical frequency to zc is found to be uj zc — 1.8T C . As can 
be seen from Eq.(Q), confinement in the radial direction 
decreases the critical temperature as well. Therefore, in 
general one would expect the existence of the critical as- 
pect ratio Ic such that the pairing is possible only in a 
trap with I > l c . 



IV. DIPOLAR BOSE GASES IN OPTICAL 
LATTICES 

In a very recent paper |34|, we have investigated the 
ground-state properties of a polarized gas of bosonic par- 
ticles (atoms or molecules) possessing a dipole moment 
and placed in an optical lattice. This system presents 
features which are novel both in the theory of quan- 
tum phase transitions, and in the context of degenerate 
quantum gases. Recently, Jaksch et al. fl2]| analyzed the 
superfluid-Mott insulator phase transition in the context 
of cold bosonic atoms with short-range interactions in 
an optical lattice. This analysis has been very recently 
confirmed experimentally JIT) . The dipole-dipole interac- 
tions, which are anisotropic and have a long range, have 
not been studied, to the best of our knowledge, in the con- 
text of the Mott insulator to superfluid phase transition. 
In addition, we shall show that the long-range character 



of the dipole-dipole interaction allows for the existence of 
not only Mott-insulating and superfluid phases, but also 
several other phases in the system. Very importantly, the 
gas of cold dipolar bosons in an optical lattice is shown to 
be a system with easily tunable interactions which may 
permit to realize all the different phases experimentally. 
The availability of such a highly controllable system may 
be crucial in answering some unresolved questions and 
controversies in the theory of quantum phase transitions 
(e.g. the existence of a yet unobserved supersolid | p5| , or 
a Bose metal at zero temperature 56 1). 



A dilute gas of bosons in a periodic potential (e.g. 
in an optical lattice) can be described with the help of 
the Bose-Hubbard model fl^] . For particles interacting 
via long-range forces the Bose-Hubbard Hamiltonian be- 
comes: 



h = j b fo + l u o n ^ - *) 



<%,J> 



1 1 

+ rurij + -U 2 2^ 



riirij + 



(8) 



<hj> 



where bi is an operator annihilating a particle at a lat- 
tice site i in a state described by the Wannier function 
w(r — Ti) of the lowest energy band, localized on this site. 
Ti is the position of the local minimum of the optical po- 
tential, and rii = b\bi is the number operator for the site i. 
In Eq. (^) only the nearest neighbor tunneling is consid- 
ered, which is described by the parameter J. The inter- 
particle interactions are characterized by the parameters 
U m , where m = \ j — i\. In particular, Uq determines the 
on-site interactions, U\ the nearest- neighbor interactions, 
U2 - the interaction between the next-nearest neighbors, 
etc. Consequently, the respective summations in Eq.(|^) 
must be carried out over appropriate pairs of sites which 
are marked by <> for the nearest neighbors, <<>> for 
the next-nearest neighbors, etc. Note that in 2D and 3D, 
i is a multi-index enumerating sites in the corresponding 
lattice. 

If only short-range interactions are present, seminal 
theoretical studies of the Bose-Hubbard model and re- 
lated theory of arrays of Josephson junctions pTj] , show 
that two different quantum phases can occur for bosons 
in a lattice, either superfluid or Mott-insulator. In the 
case of finite-range interactions additional phases as su- 
persolid and checker-board are expected [M . 

We have found the ground state of the system by 
employing a variational approach based in the so-called 
Gutzwiller Ansatz (see |l^] and references therein). We 
have shown that by modifying well-controllable pa- 
rameters a whole variety of different quantum phases 
can be achieved, including superfluid, supersolid, Mott- 
insulator, checker-board and collapse phases. This pos- 
sibility of manipulation of the corresponding quantum 
phases by modifying easily controllable external parame- 
ters (e.g. the on-site aspect ratio) makes the dipolar Bose 
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gas in an optical lattice particularly challenging for the- 
ory and experiments. In particular, possible applications 
for quantum information processing are worth mention- 
ing m this context ■ 



V. PHYSICAL REALIZATIONS OF ULTRACOLD 
DIPOLAR GASES 

One of the possible physical realizations of a gas of 
dipolar particles is provided by electrically polarized 
gases of polar molecules. This molecules can have a 
large permanent electric dipole. The creation of cold 
clouds of polar molecules has been recently demonstrated 
in experiments with buffer-gas cooling p9j, as well as 
in experiments based on deceleration and cooling of 
polar molecules by time-dependent electric fields [j60|| . 
The molecular dipole moments typically range from 0.1 
D (1 D= KT 18 charge SGS x cm) to 1 D. For bosonic 
molecules that should be sufficient for achieving BEC at 
relatively low density of the gas. For example, the dipole 
moment of the fermionic deuterated ammonia molecule 
15 ND3 is d — 1.5 D, which corresponds to an effective 
scattering length (see above, Eq.(||)) of = —1450 A. 
This is even larger than the scattering length for the 
inter-component interaction in the widely discussed case 
of fermionic 6 Li. On the other hand, for the fermionic 
molecule 14 N 16 with dipole moment d = 0.16 D, the 
corresponding scattering length — —24 A. 

Another possibility is to use magnetic atomic dipolcs 
p3|-p8[. Chromium atoms have a magnetic moment 
fi = 6^b, which is equivalent to an electric dipole mo- 
ment of d* = 6- 10~ 2 D, and an effective scattering length 
ad = —5 A. This is too small and is unlikely to result in 
any interesting fermionic effects, such as BCS pairing. 
As we mentioned in section |fi|, however, this could be 
enough to observe the influence of dipole-dipole interac- 
tions on the elementary excitations of Chromium BEC, 
provided the (still unknown) s-wave scattering length is 
not anomalously large, and it takes a value in the range 
of few tens of A. 

Permanent electric dipole moments can also be created 
by applying a high dc electric field to an atom. This 
possibility was discussed in Ref. [0, and here we only 
mention that in order to induce the dipole moment of the 
order of 0.1 D (the corresponding scattering length ~ 
-10-=- 100 A) one needs an electric field of the order of 
10 6 V/cm. Nevertheless, the influence of dipole forces on 
elementary excitation spectrum of a dipolar BEC might 
be in this case also perhaps observable using not as high 
electric fields. 

Finally, we mention the possibility of inducing a time 
averaged electric dipole moment of an atom by a stro- 
boscopic laser coupling of the ground atomic state to 
a Rydberg state in a moderate dc electric field. This 
method has been proposed in the context of bosonic gases 



in Ref. p4) . We propose to place a BEC of alkali atoms 
into a moderate static electric field. The idea is to ad- 
mix, with the help of a laser, to the atomic ground state 
the permanent dipole moment of a low-lying Rydberg 
state. Rydberg states of Hydrogen and alkali atoms ex- 
hibit a linear Stark effect in Hydrogen, for exam- 
ple, an electric field E s splits the manifold of a Ryd- 
berg states of given principal quantum number n and 
magnetic quantum number m into 2(n — \m\ — 1) Stark 
states. The outermost Stark states have (large) perma- 
nent dipole moments dp> ~ n 2 eos (with cib the Bohr ra- 
dius), and there will be an associated dipole-dipole force 
between atoms. The spacing Hlo s ~ neasE s between ad- 
jacent Stark states should greatly exceed the mean-field 
dipole-dipole interaction (and the gas temperature) in 
order to avoid interaction-induced transitions from the 
lowest sublevel to other sublevels of the manifold. 

This dipole-dipole interaction can be controlled with a 
laser []l6f . This is achieved either by admixing the per- 
manent dipole moment of the Stark states to the atomic 
ground state with an off-resonant cw laser, or by a stro- 
boscopic excitation with a sequence of laser pulses. The 
pulses should be separated by the time T, have duration 
2At -c T and area multiple of 2it. 

The field E s and the laser should be chosen such 
that they do not couple the selected lowest sublevel to 
other Rydberg (sub)states. The stroboscopic excitation 
"dresses" the atomic internal states, so that each atom 
acquires a time averaged dipole moment of the order of 
d s = n 2 eaBf, oriented in the direction of E s , where 
/ = At/T. Even though the quantity / is assumed 
to be small, the induced dipole can be rather large for 
n 1. Taking for example At = Ins, T = 10/iS, and 
n = 20, we obtain d s — 0.1D. The resulting time depen- 
dent Hamiltonian can be replaced by its time average, 
leading to Eq.(0) with d = d s . The characteristic time 
scale in Eq.([l) is provided by the inverse of the trap fre- 
quency lu^ 1 . Hence, in our case the dynamics of the sys- 
tem is described by Eq.([l]) with d — d s , if the condition 
At,T < is satisfied. 

It is important to note that, for the values of the scat- 
tering length of the order of — 10 -j- 100 A, one neverthe- 
less has a possibility to achieve the critical temperature 
of the BCS transition of the order of 100 nK at densi- 
ties n ~ 10 16 cm -3 . In the case of a single-component 
atomic Fermi gas, such densities are not unrealistic be- 
cause the inelastic processes, such as two-body colli- 
sions and three-body recombination, that usually limit 
the maximum value of the gas density, are strongly sup- 
pressed. This comes from the fact that these processes 
take place at short interatomic distances (of the order of 
tens of Angstroms) where the wave function of the rela- 
tive motion of two identical fermions vanishes due to the 
Pauli principle. 
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VI. CONCLUSIONS 

In this paper we have reviewed the properties of the 
bosonic and fermionic dipolar gases, and analyzed the 
perspectives that these systems can offer. We have 
shown, that the physics of these systems can differ quali- 
tatively in a significant way with respect to the gases in- 
teracting via Van der Waals forces. In addition, new easy 
ways of control and manipulation are possible for dipolar 
gases. For all of that, we consider that the ultracold dipo- 
lar quantum gases constitute an interesting challenge for 
theory and experiments in AMO and condensed-matter 
physics. 
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